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Abstract 

We fully elucidate the structure of the hierarchy of the con- 
nected operators that commute with the Markov matrix of the 
Totally Asymmetric Exclusion Process (TASEP). We prove for 
the connected operators a combinatorial formula that was con- 
jectured in a previous work. Our derivation is purely algebraic 
and relies on the algebra generated by the local jump operators 
involved in the TASEP. 
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1 Introduction 

The Asymmetric Simple Exclusion Process (ASEP) is a lattice model of parti- 
cles with hard core interactions. Due to its simplicity, the ASEP appears as a 
minimal model in many different contexts such as one-dimensional transport 
phenomena, molecular motors and traffic models. From a theoretical point 
of view, this model has become a paradigm in the field of non-equilibrium 
statistical mechanics; many exact results have been derived using various 
methods, such as continuous limits, Bethe Ansatz and matrix Ansatz (for re- 
views, see e.g., Spohn 1991, Derrida 1998, Schiitz 2001, Golinelli and Mallick 
2006). 

In a recent work (Golinelli and Mallick 2007), we applied the algebraic 
Bethe Ansatz technique to the Totally Asymmetric Exclusion Process (TASEP) 
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This method allowed us to construct a hierarchy of 'generalized Hamiltonians' 
that contain the Markov matrix and commute with each other. Using the 
algebraic relations satisfied by the local jump operators, we derived explicit 
formulae for the transfer matrix and the generalized Hamiltonians, generated 
from the transfer matrix. We showed that the transfer matrix can be inter- 
preted as the generator of a discrete time Markov process and we described 
the actions of the generalized Hamiltonians. These actions are non-local be- 
cause they involve non-connected bonds of the lattice. However, connected 
operators are generated by taking the logarithm of the transfer matrix. We 
conjectured for the connected operators a combinatorial formula that was 
verified for the first ten connected operators by using a symbolic calculation 
program. 

The aim of the present work is to present an analytical calculation of the 
connected operators and to prove the formula that was proposed in (Golinelli 
and Mallick 2007). This paper is a sequel of our previous work, however, in 
section [2J we briefly review the main definitions and results already obtained 
so that this work can be read in a fairly self-contained manner. In section [31 
we derive the general expression of the connected operators. 

2 Review of known results 

We first recall the dynamical rules that define the TASEP with n particles 
on a periodic 1-d ring with L sites labelled i — 1, . . . , L. The particles move 
according to the following dynamics: during the time interval [t,t + dt], a 
particle on a site i jumps with probability dt to the neighboring site i + 1, if 
this site is empty. This exclusion rule which forbids to have more than one 
particle per site, mimics a hard-core interaction between particles. Because 
the particles can jump only in one direction this process is called totally 
asymmetric. The total number n of particles is conserved. The TASEP 
being a continuous-time Markov process, its dynamics is entirely encoded in 
a 2 L x 2 L Markov matrix M, that describes the evolution of the probability 
distribution of the system at time t. The Markov matrix can be written as 



L 



1 1 




i=i 



where the local jump operator Mj affects only the sites i and i + 1 and 
represents the contribution to the dynamics of jumps from the site i to i + 1. 
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2.1 The TASEP algebra 

The local jump operators satisfy a set of algebraic equations : 

Ml = —Mi, (2) 
Mi M i+1 Mi = M i+l Mi M i+1 = 0, (3) 

[Mi,Mj\ = if \i-j\ > 1. (4) 

These relations can be obtained as a limiting form of the Temperley-Lieb 
algebra. On the ring we have periodic boundary conditions : M i+ i = Mi. 
The local jumps matrices define an algebra. Any product of the Mj's will be 
called a word. The length of a given word is the minimal number of operators 
Mi required to write it. A word, that can not be simplified further by using 
the algebraic rules above, will be called a reduced word. 

Consider any word W and call T{W) the set of indices i of the operators 
Mi that compose it (indices are enumerated without repetitions) . We remark 
that, if W is not annihilated by application of rule ([H]), the simplification 
rules (E], H]) do not alter the set Z(W), i.e., these rules do not introduce any 
new index or suppress any existing index in I(W). This crucial property is 
not valid for the algebra associated with the partially asymmetric exclusion 
process (see Golinelli and Mallick 2006). 

Using the relation ([2]) we observe that for any i and any real number 
A 7^ 1 we have 

(1 + AM;)" 1 = (1 + aMi) with a = - . (5) 

A — 1 

2.2 Simple words 

A simple word of length k is defined as a word M a mM a (2) ■ ■ ■ M a ^, where a 
is a permutation on the set {1, 2, ... , k}. The commutation rule (jlj) implies 
that only the relative position of Mj with respect to M»±i matters. A simple 
word of length k can therefore be written as W k (s2, S3, . . . , sj.) where the 
boolean variable Sj for 2 < j < k is defined as follows : Sj = if Mj is 
on the left of Mj_i and Sj = 1 if Mj is on the right of Mj_i. Equivalently, 
Wk(s2, S3, ... , Sfc) is uniquely defined by the recursion relation 

W k (s 2 , s 3 , ... , s fc _i, 1) = W k -i(s 2 , s 3 , ... , Sfc-i) M k , (6) 
W k (s 2 ,s 3 ,...,s k _i,0) = M k W / fc _i(s 2 ,s 3 ,...,s /c _i). (7) 

The set of the 2 k ~ l simple words of length k will be called For a simple 
word Wk, we define u(W k ) to be the number of inversions in W k , i.e., the 
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number of times that Mj is on the left of Mj_i : 

k 

u(W k (s 2 , s 3 , ... , s fc )) = ~ s i) • ( 8 ) 

We remark that simple words are connected, they cannot be factorized 
in two (or more) commuting words. 

2.3 Ring-ordered product 

Because of the periodic boundary conditions, products of local jump opera- 
tors must be ordered adequately. In the following we shall need to use a ring 
ordered product O () which acts on words of the type 

W = M h M i2 ...M ik with 1 < % x < i 2 < . . . < i k < L , (9) 

by changing the positions of matrices that appear in W according to the 
following rules : 

(i) If ii > 1 or i k < L, we define O (W) = W. The word W is well- 
ordered. 

(ii) If % x — 1 and i k = L, we first write W as a product of two blocks, 
W = AB, such that B = M^M^i . . . Ml is the maximal block of matrices 
with consecutive indices that contains Ml, and A = M\ .1/, , . . . .1/,, . with 
i a < b — 1, contains the remaining terms. We then define 

0(W) = (AB) = BA = M b M b+1 . . . M L M x M l2 . . . M ia . (10) 

(iii) The previous definition makes sense only for k < L. Indeed, when 
k = L, we have W = M\M 2 . . . M L and it is not possible to split W in two 
different blocks A and B. For this special case, we define 

0(M x M 2 ...M L ) = |1,1,...,1)(1,1,...,1|, (11) 

which is the projector on the 'full' configuration with all sites occupied. 

The ring-ordering O () is extended by linearity to the vector space spanned 
by words of the type described above. 

2.4 Transfer matrix and generalized Hamiltonians Hk 

The algebraic Bethe Ansatz allows to construct a one parameter commuting 
family of transfer matrices, t(X), that contains the translation operator T = 
t(l) and the Markov matrix M = t'(0). For < A < 1, the operator 
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t(X) can be interpreted as a discrete time process with non-local jumps : 
a hole located on the right of a cluster of p particles can jump a distance 
k in the backward direction, with probability A fc (l — A) for 1 < k < p, 
and with probability A p for k = p. The probability that this hole does 
not jump at all is 1 — A. This model is equivalent to the 3-D anisotropic 
percolation model of Rajesh and Dhar (1998) and to a 2-D five- vertex model. 
It is also an adaptation on a periodic lattice of the ASEP with a backward- 
ordered sequential update (Rajewsky et al. 1996, Brankov et al. 2004), 
and equivalently of an asymmetric fragmentation process (Rakos and Schiitz 
2005). 

The operator t(X) is a polynomial in A of degree L given by 

t(X) = 1 + ]T \ k H k , (12) 
k=i 

where the generalized Hamiltonians H k are non-local operators that act on the 
configuration space. [We emphasize that the notation used here is different 
from that of our previous work : t(X) was denoted by t g (X) in (Golinelli and 
Mallick 2007).] 

We have Hi — M and more generally, as shown in (Golinelli and Mallick 
2007), H k is a homogeneous sum of words of length k 

H k = ]T G(M h M i2 ...M ik ) , (13) 

l<ii<i2<...<ifc<L 

where O () represents the ring ordered product that embodies the periodicity 
and the translation-invariance constraints. 

For a system of size L with N particles only Hi, H 2 , . . . , H^ have a non- 
trivial action. Because we are interested only in the case N < L — 1 (the full 
system as no dynamics) there are at most L — 1 operators H k that have a 
non-trivial action. 

3 The connected operators Fj, 
3.1 Definition 

The generalized Hamiltonians H k and the transfer matrix t(X) have non-local 
actions and couple particles with arbitrary distances between them. Besides 
H k is a highly non-extensive quantity as it involves generically a number of 
terms of order L k . As usual, the local connected and extensive operators 
are obtained by taking the logarithm of the transfer matrix. For k > 1, the 
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connected Hamiltonians Fk are denned as 



fe=l K 



(14) 



Taking the derivative of this equation with respect to A and recalling that 
t(X) commutes with t'(X), we obtain 



Expanding t(A) -1 with respect to A, this formula allows to calculate Fk as a 
polynomial function of Hi, ... , Hk- For example F± = H\, F2 = 2H2 — Hf, 
etc... (see Golinelli and Mallick 2007). By using ffl3l) . we observe that Fk is 
a priori a linear combination of products of k local operators Mj. However 
this expression can be simplified by using the algebraic rules (0, El HI) and in 
fine, Fk will be a linear combination of reduced words of length j < k. 

Because of the ring-ordered product that appears in the expression ( TT3l) 
of the .f/fc's, it is difficult to derive an expression of Fk in terms of the local 
jump operators. An exact formula for the Fk with k < 10 was obtained 
in (Golinelli and Mallick 2007) by using a computer program and a general 
expression was conjectured for all k. In the following, the conjectured formula 
is derived and proved rigorously. 

3.2 Elimination of the ring-ordered product 

The expression A i 7 }. can be written as a linear combination of reduced 
words W. We know from formula (1131) that at most L — 1 operators Hk are 
independent in a system of size L, we shall therefore calculate Fk only for 
k < L — 1. Thus, we need to consider reduced words of length j < L — 1. 
Let W be such a word, and T(W) be the set of indices of the operators Mi 
that compose W; our aim is to find the expression of W and to calculate 
its prefactor from equation ffTol) . Because the rules ©H]) do not suppress or 
add any new index, the following property is true : if a word W appearing 
in A t(A) _1 t'(A) is such that I(W') 7^ I(W) then even after simplification, 
W will remain different from W. Therefore, the prefactor of W in X k Fk 
is the same as the prefactor of W in 



Ati(A)" 1 ^(A) where ti(A) = C n( 1 + AM i) with X(W) C X. (16) 



Because Fk commutes with the translation operator T, then for any r = 
1, . . . , L — 1, the prefactor of W — M^M^ . . . Mj . is the same as the prefactor 



00 



^A fe F fc = At(A)- 1 f(A). 



(15) 



fc=i 
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of T r MT~ r = M r+il M r+i2 . . . M r+ij . Furthermore, any word W of size k < 
L — 1 is equivalent, by a translation, to a word that contains Mi and not 
Ml ■ indeed, there exists at least one index i such that io £ Z(W) and 
(zo + 1) G 1{W) and it is thus sufficient to translate W by r = L — i . 

In conclusion, it is enough to study in expression (TT3Ti . the reduced words 
W with set of indices included in 

J* = {1,2,...,L-1}. (17) 

Because the index L does not appear in X*, the ring-ordered product has a 
trivial action in equation (1161) and we have 

t T *(\) = (1 + AMi)(l + AM 2 ) . . . (1 + AM L _0 . (18) 

We have thus been able to eliminate the ring-ordered product. 

3.3 Explicit formula for the connected operators 

In equation (Tig]) , differentiating tj*(A) with respect to A, we have 

L-1 

t' T * (A) = £(l + AM 1 )...(l + AM J _ 1 )M,(l + AM m )...(l + AM L _ 1 ). (19) 
i=i 

Using equation (0) we obtain 

£z*(A) -1 = (l + aMi_ 1 )(l + aM i _ 2 )...(l + aMi), with a = -A- . (20) 

A — 1 

Noticing that A(l + aMj)Mj = we deduce 

A MA)" 1 4(A) = (21) 

L-1 

-a ^(1 + «M L _0 . . . (1 + aMi+OM^l + AM m ) . . . (1 + AM L _ X ) . 

i=l 

The ith term in this sum contains words with indices between i and L — 1. 
Because we are looking for the words that contain the operator M\, we must 
consider only the first term in this sum, which we note by Q 

Q = - a (l + aM w ) • • • (1 + aM 2 )M 1 (l + AM 2 ) . . . (1 + AM^) . (22) 

In the appendix, we show that 

Q = R 1 + R 2 + ... + R L ^ , (23) 
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where Ri is defined by the recursion : 

R x = -aM 1 , (24) 
Ri = \Ri-xMi + aMiRi-! for % > 2 . (25) 

To summarize, all the words in Y^kLi ^ k Fk that contain M\ and not Ml are 
given by Q = Ri+R 2 + - ■ .+Rl-i- From the recursion relation (|25|) we deduce 
that Ri is a linear combination of the 2 l_1 simple words Wi(s2, S3, Sj) 
defined in section 12.11 Furthermore, we observe from f)25p that a factor A 
appears if Si = 1 and a factor a = A/(A — 1) appears if s, = 0. Therefore, 
the coefficient /(W 7 ) of = Wi(s2, S3, ... , Sj) in Q is given by 

\i 00 / \ 

/ w = (-ir (1 _ A)M+1 = E (7 j Ai+j (26) 

where z is the length of W and u = u{W) is its inversion number, defined in 
equation (jHJ). We have thus shown that 

Q = E E /W^ = E E ^E(-ir w f uW+ 'V^ ; (27) 

i=i iyGW t i=i vveWi j=o V J / 

where Wi is the set of simple words of length i. 

Finally, we recall that the coefficient in J2kLi ^ k Fk of a reduced word W 
that contains Mi and not M^ is the same as its coefficient in Q. Extracting 
the term of order A fe in equation (1271) we deduce that any word W in Fj, that 
contains M\ and not Ml is a simple word of length i < k and its prefactor 
is given by 

The full expression of Ffc is obtained by applying the translation operator 
to the expression fl27|) ; indeed any word in F k can be uniquely obtained by 
translating a simple word in F k that contains Mi and not Ml- We conclude 
that for k < L, 

F k = T±J2 (-ir^( k - t+u(W) )w , (28) 
i=i weWi V k 1 ) 

where T is the translation-symmetrizator that acts on any operator A as 
follows : TA = J2i=o T l A T~ l . The presence of T in equation ( J28l) insures 
that Ffc is invariant by translation on the periodic system of size L. All simple 
words being connected, we finally remark that formula (1281) implies that Fk 
is a connected operator. 
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4 Conclusion 

By using the algebraic properties of the TASEP algebra (EJHJ), we have derived 
an exact combinatorial expression for the family of connected operators that 
commute with the Markov matrix. This calculation allows to fully elucidate 
the hierarchical structure obtained from the Algebraic Bethe Ansatz. It 
would be of a great interest to extend our result to the partially asymmetric 
exclusion process (PASEP), in which a particle can make forward jumps 
with probability p and backward jumps with probability q. In particular, we 
recall that the symmetric exclusion process is equivalent to the Heisenberg 
spin chain : in this case the connected operators have been calculated only 
for the lowest orders (Fabricius et al., 1990). This is a challenging and 
difficult problem. In our derivation we used a fundamental property of the 
TASEP algebra : the rules (EJHj) when applied to a word W either cancel 
W or conserve the set of indices X(W). The algebra associated with PASEP 
violates this crucial property because there we have Mj M i+1 Mj = pq Mj. 
Therefore the method followed here does not have a straightforward extension 
to the PASEP case. 



Appendix: Proof of equation (1231) 

Let us define the following series 

Qi = -aMi , (29) 
Qi = (1 + aMi)Qi_i(l + AM») for i > 2 . (30) 



We remark that Q defined in equation (12211 is given by Q = Ql-i- Let us 
consider Ri defined by the recursion (1231) . The indices that appear in the 
words of Qi and Ri belong to {1,2, ... Therefore, we have 

[R j ,M i ]=0 for j<i-2, (31) 

because the operators Mi, M 2 , . . . , Mj that compose Rj commute with Mj. 
From equations (13T1) and (3J, we obtain 

(1 + aMi)Rj(l + AM,) = Rj hrj<i-2. (32) 

Furthermore, from (123|) . we obtain 

MiRt-iMi = A .\ /,/.', ,.U, : M, + a.\l,M, .R, ,.U, . (33) 

Because Mj commutes with R4-2, we can use the relation MiMi_\Mi = to 
deduce that 

M i Ri_ 1 M i = 0. (34) 
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Using equation f[3"4l . we find 

(1 + aM^Ri^l + AM;) = + XRi-iMi + aMiRi_i = + R { . (35) 

From equations (1321 and fl35|) . we prove that the (unique) solution of the 
recursion relation flHUl) is given by equation f[2"3"j) . Q; = i^i + i? 2 + • • • + Ri- 
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